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, We study the problem of motion of a relativistic, ideal elastic solid with 

qh| free surface boundary by casting the equations in material form ("La- 

grangian coordinates"). By applying a basic theorem due to Koch, we 
bj). prove short-time existence and uniqueness for solutions close to a trivial 

solution. This trivial, or natural, solution corresponds to a stress-free body 
in rigid motion. 



Abstract 



1 Introduction 

In the problem of motion for classical continua with free surface boundary, despite 
its obvious physical relevance, there is a surprising scarcity of rigorous results. 
A recent review of known results with and without gravity can be found in [14]. 
In the case of a nonrelativistic, non-gravitating perfect fluid local wellposedness 
has only been proved very recently by Lindblad [12]. In the present paper we 
consider the analogous problem for relativistic elastic solids without self-gravity. 
The nonrelativistic case of our results (which is in principle known - see [15] - 
although even there a detailed treatment seems to be lacking) can be proved by 
similar methods. The nonrelativistic case for incompressible materials has, by 
different methods, been treated in [7]. 

The plan of this paper is as follows. In Sect. 2 we describe elasticity theory on 
an arbitrary curved background as a Lagrangian field theory. Sect. 3 is devoted 
to the notion of "natural state", that-is-to-say a solution to the elastic equations 
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corresponding to a configuration with zero stress. The existence of such a solution 
requires the elastic body to move along a geodesic, timelike Killing vector. We 
require the Killing vector to be also hypersurface-orthogonal (thus forcing space- 
time to be "ultrastatic")0. In the case where the metric on the space orthogonal 
to the Killing vector is flat, we are in Special Relativity. In Sect. 4 we perform a 
3+1 decomposition of the elastic equations corresponding to the natural space- 
time splitting afforded by the Killing vector. We then write the elastic equations 
in "material" form (often called "Lagrangian representation"). This means that 
the maps making up the elastic configuration space - which go from spacetime 
to the 3 dimensional "material space" - are replaced by time dependent maps 
("deformations") from material space into physical space. In the material rep- 
resentation the boundary of the body is fixed, namely the boundary of material 
space. In Sect. 5 we state a corollary to the theorem of Koch [11], which is the ver- 
sion of existence theorem we are using. In Sect. 6 we state our basic constitutive 
assumption. This assumption, which is satisfied by elastic materials occurring in 
practice, implies the validity of the conditions in the Koch theorem. From this 
one concludes the main result, which is stated in Sect. 6. In the appendix we show 
that the corner conditions on the boundary, which initial data have to satisfy in 
order for the time evolved solution to be a classical solution, can be satisfied for 
a large class of initial data. 

We should emphasize that our requirement of the existence of a global relaxed 
state, and the associated restriction on the background spacetime would be un- 
necessary if we were dealing with well-posedness in the absence of a boundary, or 
inside a region of the body causally disconnected from the boundary. However 
in the presence of the boundary we have to ensure corner conditions using an 
implicit function argument, and for this we use the existence of a state - the 
global relaxed one - which satifies these corner conditions to all orders. In this 
sense our strong restriction on the existence of a global relaxed configuration is 
perhaps not completely essential, but very convenient. 

We can state the central result of this paper as follows: For initial data close to 
initial data for the natural deformation in the appropriate function space, there 
exists, for sufficiently short times, a unique solution to the elastic equations. This 
solution depends continuously on initial data, in particular tends to the natural 
deformation, when the initial data tend to that of the natural deformation. 

2 The theory 

We treat elasticity as a Lagrangian field theory in the manner of [2] or [4], see 
also [10]. In the language of standard elasticity this means that the material in 
question is "hyperelastic". The dynamical objects of the theory are furnished by 

""^In Special Relativity a geodesic, timelike Killing motion has to be a time translation, so 
the above assumption is superfluous in that case. 



2 



sufficiently regular maps / sending a closed region S = S U dS of spacetime M. 
onto B = B U dB, with B a bounded domain in with smooth boundary dB 
called the "body" or material space. The domain B is to be thought of as an 
abstract collection of points ("labels") making up the elastic continuum, and / 
is to be thought of as the back-to-labels map sending spacetime events to the 
particles by which they are occupied. We endow B with a volume form Q which 
is smooth on B. Coordinates on Ai are written as x^, with /i, ly etc ranging from 
to 3 and coordinates on B as with A, B etc going from 1 to 3. The metric 
g^^i, is also taken to be smooth. There are additional requirements on the maps 
/. Namely the equation 

= 0, g,,u^v!' = -1 (1) 

should have a solution which is unique up to sign. Thus particles move along 
trajectories in with unit tangent u^. Furthermore the function n on S defined 
by the equation 

nABcU{^))f^A^)f''A^)fA^) = n{x)e,,^p{x)uP{x) (2) 

should be everywhere positive on S. It is well known that ([2]) implies that the 
continuity equation 

V^inu^') = (3) 

is identically satisfied. The equations of motion are the Euler-Lagrange equations 
of the action 

S[f,df]= J^p{f,df,x){-gy^ d'x. (4) 

The function p is required to be smooth in all its arguments. It plays the double 
role of being the Lagrangian as well as energy (i.e. rest mass plus internal elastic 
energy) density of the material. Diffeomorphism covariance demands that p be 
of the form 

p{f,df,x) = p{f,H^^)\HAB=hAB^Qf^^), (5) 

where h^^{df,x) = f^^^{x)f^ ^y{x)g^^{x). Thus p depends only on points X on 
B and positive definite contravariant tensors H"^^ (see [2]). (In nonrelativistic 
elasticity materials governed by such a Lagrangian are said to satisfy the condition 
of "material frame indifference".) By abuse of notation we henceforth omit the 
hat from p. It is useful to know that the Euler-Lagrange equations are equivalent 
to divergence-lessness of the symmetric energy momentum tensor. More precisely, 
there holds the identity (see [2]) 

- V.T/ = £^/^,^, (6) 
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where 

^/^- = 2^ - P9^.u (7) 

and 

The relation ([6]) shows that of the four conservation laws Vj^T^'^ = only three are 
independent. This, of course, is due to the fact that within the present formalism 
the continuity equation is an identity. 

We now turn to boundary conditions. The boundary conditions usually consid- 
ered in nonrelativistic elasticity are either the so-called "boundary conditions of 
place" - where, in our language, the set f~^{dB) is prescribed, or "traction bound- 
ary conditions" - where the normal traction, i.e. the components of the stress 
tensor normal to this surface are prescribed. When the normal traction is zero, 
one speaks of "natural" boundary conditions: these are the boundary conditions 
appropriate for a freely floating elastic body considered in the present work. They 
are employed e.g. in geophysics for the elastic motion corresponding to seismic 
waves, the free boundary in question corresponding to the surface of the earth 
(see [1] [§). While one could in the present framework in principle consider all the 
above boundary conditions, it is interesting to observe that these conditions - 
except for the natural condition - become inconsistent once one couples to the 
Einstein equations. Namely the standard junction conditions, together with the 
Einstein equations, imply that T^'^n^\f-i(^g^) be zero, where is the conormal 
of f~^{dB). However these conditions are precisely the natural boundary condi- 
tions. To see this, we have to first compute the right hand side of Eq.(l7l). The 
result is 

Tf,u = pu^u^ + t^,y, (9) 
where tfj,^, is the (negative) Cauchy stress tensor given by 

V = nTABf^,j'',. (10) 

and we have written p as 

p = ne (11) 

and Tab = 2 qI^ab ■ Note this makes sense, since n - whence e - is a function 
purely of / and H^^, as is apparent from the identity 

6n' = H^^'H^^'H^^'Qabc^A'B'c. (12) 

The function e is the relativistic version of the " stored-energy-function" of stan- 
dard elasticity. The quantity tab corresponds to the negative of the "second 



^Needless to say, these studies are confined to the linear approximation in which the problem 
of the free boundary disappears. 
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Piola-Kirchhoff stress" of nonrelativistic elasticity. Now to the boundary condi- 
tions, namely 

T^j,'"n^\f-i(dB) = 0. (13) 

It follows from the tangency of to the inverse images of points of B under /, 
that (tt,n) = 0. Consequently the contraction of Eq.( fT3i) with u^, using IQ, is 
identically satisfied. The remaining components yield 

rABf'',^9'"'n,\f-i^9B) = (14) 

Equation (IMI) will turn out to be a Neumann-type boundary condition, but it 
has a free (i.e. determined-by-/) boundary. 



3 Natural configuration 

We assume there exists a contravariant metric Hq^{X), smooth and positive 
definite on B, such that 



de 

\(x,HAB^H^B(^x)) = eo(^) > 0, 



= (15) 

on B. The quantity H^^ measures local distances in the relaxed body. The 
function eo(-^) is the density of rest mass in the stressfree configuration. 
If there exists a configuration /q such that 

H^^'iUx)) = h^''{dMx),x) = ft,{x)f^A^)r{x) =: hf{x), (16) 

this map /o is called a natural or relaxed configuration. Note that spacetime has 
to be special in order for a natural configuration to exist. Namely it follows from 
( IT6|) that h^^ has to be constant along «q, the four velocity associated with /o, 
from which one infers that Uq is a Born rigid motion (see [16]), i.e. 

^uo{9fiu + Uo^Uqi,) = (17) 

One deduces from fl6H71l9l) and f|T5l) that /o solves the equations of motion 8,a = 
if and only if this rigid motion is geodesic. This in turn implies that Uq is Killing. 
In this work we will assume that Uq is in addition irrotational. Using coordinates 
in which u^d^ = dt and u^^dx^ = —dt, the spacetime metric has thus to be of 
the form (i, j = 1, 2, 3) 

g^^dx'^dx" = -dt^ + gij{x^)dx''dxK (18) 

(In particular: when the positive definite three metric gij is fiat, we are in 
Minkowski spacetime.) With this choice of coordinates, the natural configu- 
ration /o is of the form fo{t,x) = fo{x) with fo{x) a smooth, invertible function 
S O {t = 0} ^ B where, using coordinates on M and B in which eijk and ^abc 
are both positive, det{fQ^i) is positive on 5 fl {t = 0}, the latter condition guar- 
anteeing that Uq is positive and is positive definite on 5 fl {t = 0}. The 
inverse of fo{x) will be denoted by $o(X). 
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4 3 + 1 and the material picture 

In this section we switch to the material description of the elastic medium by 
changing to "Lagrangian" coordinates. In a first step we replace the four velocity 
corresponding to the configuration / by the coordinate velocity f ^(9^ = dt + 
V'di, which is a multiple of u^. Using ([1]) we find that V" is given in terms of f"^ 
by (note that f^^i is non-degenerate by construction) 

where a dot denotes partial differentiation w.r. to t. It follows from fll8p that 

/i^"" = f'',,f'',u9'"' = f^,f, {g'' - V'V). (20) 

The timelike nature of (equivalently: the positive definiteness of h^^^ requires 
that 

\V\^ = g^,Y'V' <\ (21) 

We also see that 

n = A; (1 - 1^12)1/2, (22) 

where k is defined by 

k eijk = f^,if^ ,jf'^ ,k ^ABC (23) 

with eijk being the volume element of gij. Consequently the action (jll) takes the 
form 

S= [{1- \V\^y/\k{det Qijfl^dtd'x (24) 



We now pass over to the material representation. By this we mean that configu- 
rations are replaced by "deformations" = $*(f, X) defined by 

f\tMt.m=X^ (25) 
with det{(^\A) positive in B. It follows from fl25l) that 

f^,{tMt.X))^\B{t,X)=5^B. ^\A{t.X)f^,{tMt.X)) =6', (26) 

and 

f^it, $(t, X)) + f^,{t, ^t, X))^'{t, X) = (27) 
and, from (^3) and that 

V'{t,^{t,X)) = ¥{t,X) (28) 

In order to study the field equations in the material representation it will be 
extremely convenient to change representation directly in the action fl2^ . Using 
(|23|) there holds 

S = [ e-f-^dtd^n (29) 

J{t}xB 



with 



7-1 = (1- |V|2)2 (30) 

and d^n = nABciX) dX"^ A dX^ A dX'^ = d^X. Now the field equations take 
the form 



a 



dL 



dL 



dL 



inS 



(31) 



where L = e7~ fia for functions $ such that det{^\) is positive in B. Note that 
fl depends only on X^, gij depends only on and 7"^ depends only on ($\ $■'). 
Finally h^^ , i.e. 



(32) 



where F^i{(^\B) is the inverse of $\a, depends on ($*, $'^,a), so that e depends 
on (X^, $\ •i)-', ^'^^s). Using these facts, together with the identity 



dh 



AB 



we find that 



and 



dL 



dL 
9$S 



de 



QHBC 



HDE^hDE 



(—) 

\dH^^ J 



B 



HCD=hCD 



(33) 



(34) 



(35) 



From (IMI) we infer that the boundary conditions (fT4|) are simply equivalent to 

dL 



d^' 



0, 



(36) 



{t}xdB 



where ua is a conormal of dB. 

We now turn to the natural deformation, which is the material version of the 
natural configuration described in Sect. 3. In the chosen foliation it is time- 
independent, namely of the form 



^\t,X) = %{X), 
where $o(/o(3;)) = Furthermore we have that 

F^($o',c)F^($o',d)^7^^($oW) = i/o^^(X), 



(37) 



(38) 



in short: F^^i Fq^ j qq''^ = Hq^ . We know from Sect. 3 that the field equations 
and the boundary conditions are identically satisfied in the stress-free state. In 
the present context this fact takes the form 



dL 



0, 



dL 
7 



0, 



dL 



0, 



(39) 



where the subscript means evaluation in the stress- free state $*) = ($o*) 0). 
To derive Eq.([39D we have used (fT5ll3ll35D . 

The causahty properties of this system are essentially governed by the nature of 

the coefficients M"jj of and M^^ij of ^\ab- In order to be able to apply 
the theorem of Koch [11], we will need a "negativity" condition on M**jj and a 
certain positivity (" coerciveness" ) condition on M^^ij. We easily find that 

Kij = T^T^ = ^0 90^J (40) 



For the coefficient of ^ab, using (!34l) . observe that 



where 



^OABCD 



( ) 



(42) 



5 Koch theorem 



For convenience we state here a corollary of the theorem in [11], which is the 
precise statement we need. 

Let $*(t,X) be maps from {t} x B to M^, where B is an open set in M.^ with 
smooth boundary dB. We are given a system of the form of 

da3"^i = Wi in {t} X B, (43) 

where a runs from to 3 and with 3^°'i and Wj all being smooth functions of 
{X^, $*, ,a) on ;B X X R^^, together with the boundary conditions 

^^"inalitjxdB = ^"^inAlltjxdB = 0. (44) 

We make the following further assumptions: 

(i) Symmetry: there hold the symmetries M^^ij = M^^'ji where M°'^ij = 

(ii) Static solution: There is given a time independent function <l>Q(f, X) = 
%{X) G C^{B) satisfying 

3'^, = 0, Wo, = (45) 

and the following estimates 

(iii) Time components: 

M^\jr]'r]^ < -K|r/P in B (46) 
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for a positive constant 
(iv) Space components: 

/ Mo%- 6^\AS^\Bd'X+\mlHB) > ^imilHB) (47) 
Jb 

for cr > and all 5$ G C°°{B), where i/** denotes the Sobolev space if'*'^. We re- 
mark that condition (H71) implies the " strong ellipticity" or " Legendre-Hadamard" 
condition, namely that MQ^^jNANsk^k^ > in i3 for non-zero N^^ k\ This latter 
condition is the relevant one for wellposedness in the pure initial value problem 
(see [9]). 

Theorem: Let ($(0), $(0)) lie in a small neighborhood of ($o,0) in H''^^{B) x 
H^{B),s > 3 and satisfy the corner conditions of order s in the sense that 
dl3'fnA\{t=o}xdB is in H'-''{B) n Hl{B) for < r < s |. Then there exists, 
for sufficiently small to; a unique classical solution $ of fl43|44l) in C^([0,to) x 
with ($(0), <i>(0)) as initial data. Furthermore G L'^{B) for < r < s + 1. 

In the last expression 9'"$ denotes all partial derivatives of order r. The evolved 
solution stays close to the static solution $o in the sense that ^if)) remains 

close to ($0,0) in H'+\B) x H'{B) and 9^$(t) remains close to 8''% G L'^{B) 
for < r < s + 1. (This last statement is the "stability" part of the theorem.) 

6 Main result 

We now add our final constitutive assumption, sometimes called "pointwise sta- 
bility" see [13], namely that 

LoABCD iV^^iV^^ > a {HocaHobd + HocbH^ad) iV^^iV^^ in B (48) 

where cr is a positive constant which only depends on the choice of coordinates 
and with Hqab being the inverse of Hq^. Condition fHS]) is a convexity condition 
on /9 as a function of H^^. We note that this in general does not imply convexity 
with respect to the deformation gradient (see [5]). An important special 
case is where 

4:€oLoABCD = ^HqabHqcd + ^{HqcaHqbd + HqcbHqad) (49) 

and 

^(X)>0, 3\{X) + 2fi{X) > inB (50) 

The quantities /i and A, when they are independent of X, are the Lame constants 
of homogeneous, isotropic materials. We now invoke the Korn inequality (see 

•^This condition says that {dfS^'^i) nA\{t=o}xdB is zero for < r < s In the appendix we 
show that it can be met by a suitable choice of odd-order normal derivatives of ($(0), $(0)) on 
dB. 
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e.g. [8] or [6]) of which we need a shght generahzation due to [3], in the following 
form: Let "if^ be vector field on {B,Hqab) in some chart. Let Li^^) be defined 
by L{'if)AB = 2-ffoc(A^*",B), i-e. the Lie derivative of Hqab w.r. to \1/ to leading 
derivative-order in Then there is a positive constant a' such that 

\\Lm\l^iB) + \\nUB)>^'\\nmm- (^i) 

We now assume coordinates X m. B to be chosen so that $o is the identity map. 
Combining (14T|) with (148|) and using (I5T1) . there follows condition (iv), i.e. that 

/ M^'',^5^\A5^\Bd'X + \m\l^^B) > (^"\m\m(B) (52) 
Jb 

for positive a". The validity of the condition (iii) in the Koch theorem is immedi- 
ate from Eq. fl40l) . Furthermore the validity of (ii) has been checked in Sect. 4. The 
validity of (i) is obvious from the variational character of the equations. Thus all 
assumptions of the Koch theorem are met. 

It remains to check the determinant condition det{^^^A) > in i3. But, when 
($(0), $(0)) is close to ($o, 0) in H'+\B)x H'{B), s > 3, this immediately follows 
by Sobolev embedding and the positivity of (iet($Q,A) in B. Now the precise form 
of our final statement can be read off from the Koch theorem in Sect. 5. 
Stated somewhat informally, our final result is as follows: 

Theorem: Let there be given a volume form Q{X^), a spacetime metric g^u{x^) 
of the form (fT8|) and an internal energy e{X^, H^'^) satisfying (fT5|) . all smooth 
functions of their respective arguments. Suppose there exists a smooth natural 
(i.e. stress-less) configuration such that the corresponding vector field Uq is the 
static Killing field dt- Also suppose that the elasticity tensor for this natural 
configuration satisfies the pointwise stability condition psjl . Let the initial data 
($(0), $(0)) for the deformation be close to those for the natural deformation and 
satisfy the corner condition of the appropriate order. Then there exists, for these 
initial data, a solution $(t) of the dynamical equations flHTl) with boundary condi- 
tions (!36|) for small times. This solution remains close to the natural deformation. 

We end with a remark on the possible generalization of the results presented 
here to the case of an elastic body (or bodies, if there are several) which are 
self-gravitating. The presence in G.R. of constraints and gauge freedom will then 
make things more complicated. Furthermore the fact that the elastic equations 
are written in material form, but the Einstein equations are equations on space- 
time, means that one is now not dealing with a system of partial differential 
equations in the strict sense. 

Acknowledgement: We are grateful to Bernd Schmidt for very useful remarks 
on the manuscript. One of us (R.B.) thanks Tom Sideris for pointing out the 
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7 Appendix: corner conditions 

We study corner conditions for the system (H3|l with the boundary conditions 
dBj). Suppose that ($(0) lag, $(0) lag) is sufficiently close to ($o|aB,0). In order 
for obtaining a solution of the equations of motion, one has to be able to satisfy 
conditions on the behaviour of certain normal derivatives of ($(0), $(0)) on dB, 
which guarantee that, not only S^^iUAldB = at f = 0, but also a sufficient 
number of time derivatives at t = of that condition is satisfied. (It is understood 
that higher-than-second-order time derivatives are eliminated in terms of spatial 
and lower-order time derivatives, using the equations of motion: this by virtue of 
the negativity of M^ij is always possible for ($(0)|aB, ^(O)!^^) close to {^o\dB, 0).) 
First recall that ^Q{t,X) = ^q{X) identically satisfies the equations of motion 
and the boundary condition. Now we check that the corner condition of order 0, 
i.e. the undifferentiated-in-time Eq.(!36l) can be satisfied by a choice of (9„$(0)|9b 
(where 9„ denotes any derivative transversal to dB). To see this notice first 
that, by the above observation, {^o\dBi^) solves the order-0 corner condition. 
Furthermore 

Then the result follows from the (finite-dimensional) implicit function theorem 
using that, by virtue of (l48l) . the quadratic form 

M^^i, uaUbU (54) 

is nonsingular, which in turn follows from the strong ellipticity condition (see 
Sect. 5). One now performs the process of taking (s > 1 say) time derivatives of 
the boundary condition and eliminating 9["$(0)|ag for s + 1 > m > 1. Using 
that is zero for ($(0)1^^, $(0)|aB) = {^o\dB, 0), the s-th order corner condition 
becomes an equation of the form 

M^^ij uaUb <9^$-'' (0) \aB = lower order (55) 

for odd s and 

M^^ij hatib (0) 1 9B = lower order (56) 

for even s, where " = lower order" means expressions which depend on lower-order, 
even-numbered normal derivatives of ($(O)|a0, $(0)|ag), modulo terms which de- 
pend on normal derivatives of the same order, but which are zero when 
($(0)|aB, $(0)19^) = ($o|9H,0). The equations fl55|56p are identically satisfied if 



-riATlB 



(53) 



aB 
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('^'(0) las, $(0) las) and their normal derivatives appearing on both sides are re- 
placed by those of ($o|aB;0). Using the fact that M^^ij nAUB depends only on 
first derivatives of we are now able to recursively solve the corner conditions, 
with $(0)|aB and $(0)|ae and their even normal derivatives given arbitrarily, pro- 
vided that ($(0)|oe, $(0)|9e) is close to ($o|oB) 0). We have to - and by the above 
can - choose a large class of initial data ($(0), <i'(0)) close to ($0i0) so that the 
corner conditions are fulfilled for arbitrary order and the negativity of M**jj and 
the coerciveness of M^^j^ are satisfied. 
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